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ABSTRACT 

Wc present a Hubble Space Telescope/Wide Field Planetary Camera 2 weak-lensing study of A520, 
where a previous analysis of ground-based data suggested the presence of a dark mass concentration. 
We map the complex mass structure in much greater detail leveraging more than a factor of three 
increase in the number density of source galaxies available for lensing analysis. The "dark core" that is 
coincident with the X-ray gas peak, but not with any stellar luminosity peak is now detected with more 
than 10 (7 significance. The ~1.5 Mpc filamentary structure elongated in the NE-SW direction is also 
clearly visible. Taken at face value, the comparison among the centroids of dark matter, iiitracluster 
medium, and galaxy luminosity is at odds with what has been observed in other merging clusters with 
a similar geometric configuration. To date, the most remarkable counter-example might be the Bullet 
Cluster, which shows a distinct bow-shock feature as in A520, but no significant weak-lensing mass 
concentration around the X-ray gas. With the most up-to-date data, we consider several possible 
explanations that might lead to the detection of this peculiar feature in A520. However, we conclude 
that none of these scenarios can be singled out yet as the definite explanation for this puzzle. 
Subject headings: cosmology: observations — dark matter — galaxies: clusters: individual (A520) — 
galaxies: high-redshift — gravitational lensing: weak — X-rays: galaxies: clusters 



1. INTRODUCTION 

Numerical simulations have successfully demonstrated 
that galaxy clusters form at the intersections of filaments 
and they grow by accreting other clusters/groups pre- 
dominantly along the filaments. Because of the dom- 
inance of this virtually one-dimensional accretion, the 
cores of the galaxy clusters are subject to frequent near 
head-on collisions and thus are dynamically active. 

One of the key tools for studying merging clusters is 
the comparison among the distributions of the three clus- 
ter constituents: galaxies, hot plasma, and dark matter. 
For example, in merging clusters the intrachister mcdiiim 
suffers from ram pressure and lags behind galaxies and 
dark matter (e.g., Clowe et al. 2006; Jee et al. 2005a, 
2005b), which are believed to be effectively coUisionless. 
The contrast between collisional and coUisionless compo- 
nents becomes highest when we observe merging clusters 
at their core pass-through, when both the medium veloc- 
ity and the effect of ram pressure stripping are largest. 

The Bullet Cluster (Clowe et al. 2006) provides a re- 
markable example of the separation of cluster compo- 
nents. It possesses textbook examples of both bow shock 
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and clear offsets between dark matter/galaxy and X-ray 
numerical simulation of the cluster (Springel & 
Farrar 2007) suggests that the subcluster (Bullet) is mov- 
ing away from the main cluster at ^ 2700 km s^^ with 
respect to the main cluster. As the mass reconstruction 
is blind to the distribution of the cluster galaxies, the 
agreement of the mass clumps with the cluster galaxies 
and the offset of the X-ray gas from these arc strong ev- 
idence for coUisionless dark matter (Clowe et al. 2006). 

Another merging cluster showing a comparably re- 
markable bow shock feature is A520 at z = 0.201 (Marke- 
vitch et al. 2005). However, a weak-lensing study by 
Mahdavi et al. (2007, hereafter M07) finds a very per- 
plexing mass structure. In addition to three mass clumps 
that are spatially coincident with the cluster galaxies, the 
mass reconstruction also shows a significant (~ 4cr) mass 
peak that lies on top of the X-ray luminosity peak, but 
that is largely devoid of bright cluster galaxies. Such a 
peculiar substructm-e, referred to as a "Dark Core" in 
M07, does not appear in the weak-lensing mass map of 
the Bullet Cluster, which is believed to be at a similar 
merging stage. M07 discussed several possible explana- 
tions such as background cluster contamination, bright 
galaxy ejection, line-of-sight (LOS) structure, or vio- 
lation of the upper limit. One extreme interpretation 
would be that the dark matter particle's collisional cross 
section might be considerably larger than the upper limit 
{<Jm/mDM < f cm^ g~^) derived from the Bullet Cluster 
(Markevitch et al. 2004). Recently, independent results 
on the dark matter cross section supporting the collision- 
less nature have been reported from the studies of A2744 
(Merten et al. 2011), MACS J0025.4 1222 (Bradac et al. 
2008), DLSCL .10916.2+2951 (Dawson et al. 2011), halo 
ellipticities (Feng et al. 2010; see also Miralda-Escude 
2002), etc. However, WiUiams & Saha (2011) claim sig- 
nificant detection of light/mass offsets in A3827, which 
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can be interpreted as evidence for coUisional dark matter. 

In this paper, we present our weak-lensing study of 
A520 with Hubble Space Telescope (HST) / Wide Field 
Planetary Camera 2 (WFPC2) images. This is a critical 
follow-up study because high resolution imaging provides 
more usable galaxies for lensing analysis, and thus en- 
hances both the significance and resolution of the mass 
reconstruction. These data verify the reality of the dark 
core and also refine the mass estimates of the substruc- 
tures, which provide key ingredients for detailed numer- 
ical simulations. 

We assume {^lM^^h,h) = (0.3,0.7,0.7) for cosmology 
unless explicitly stated otherwise. This gives a plate scale 
of -3.3 kpc/" at the redshift (z = 0.201) of A520. AU 
the quoted uncertainties are at the la (~68%) level. 

2. OBSERVATIONS AND DATA REDUCTIONS 

The cluster A520 was observed in the F814W filter 
with WFPC2 in nine contiguous pointings covering ap- 
proximately the 8f5 X 5.'5 NE-SW elongated mass struc- 
ture. Each pointing was dithered four times with a total 
integration of 4, 400 s. Our data reduction starts with 
the products of the STScI CALWP2 (Gonzaga et al. 2010) 
task. First, we "drizzled" (Fruchter & Hook 2002) each 
exposure individually to obtain rectified coordinates of 
objects to be used for alignment. Then, the alignment 
information was fed into the MultiDrizzle (Koekemoer 
et al. 2002) software to remove cosmic rays and create a 
distortion-free mosaic image. We chose a Gaussian inter- 
polation kernel with a final pixel size of O'.'OG in order to 
reduce otherwise apparent aliasing due to the undersam- 
pled point-spread function (PSF) . We discarded the data 
on the Planetary Camera (PC) for simplicity in data re- 
duction, and uniformity in angular resolution during the 
subsequent analysis. 

The charge transfer inefficiency (CTI) of WFPC2 is 
potentially important because the deferred charge release 
smears and elongates shapes of objects, which must be 
distinguished from the distortion by gravitational lens- 
ing. We measured the CTI by examining the elliptic- 
ity of sub-seeing features (cosmic rays, hot/warm pixels, 
etc.), which are non-astronomical objects and therefore 
not subject to the PSF of the instrument (Jee et al. 2011; 
hereafter JH)- We refer readers to Jll for technical de- 
tails. Figure [lla) demonstrates how a hypothetical cir- 
cular WFPC2 PSF placed at different positions would 
be distorted by CTI. Both parallel and serial CTI effects 
are found to be important while it is apparent that the 
effect in the parallel direction is much stronger. 

The impact of this uncorrected CTI on our mass re- 
construction is shown in Figures [T|b) andjljc). Since no 
gravitational lensing signal is present here, the mass re- 
construction would indicate potential impact of WFPC2 
CTI on our mass reconstruction if no CTI correction were 
made. The periodic variation due to the observational 
dither pattern is clear. However, even the maximum am- 
plitude of this uncorrected systematics for the average 
source galaxy is an order of magnitude smaller than the 
average lensing signal. Therefore, the substructures in 
our mass reconstruction cannot be affected by any im- 
perfect CTI correction. 

We construct a WFPC2 PSF library from archival im- 
ages containing dense stellar fields (Jee et al. 2007a). 
Object ellipticity is determined by modeling the object 



surface profile with an elliptical Gaussian in the absence 
of the telescope seeing and the CTI. This is implemented 
by convolving the elliptical Gaussian with the model 
PSF (tweaked to simulate the CTI effect) prior to fitting 
(Jll). We select our source galaxies whose F814W mag 
is fainter than 22 while removing s puri ous objects and 
cluster galaxies defined in M07 (see \ 3.3 for details). The 
resulting number density is ~93 galaxies arcmin"^, ap- 
proximately a factor of 3-4 higher than the number den- 
sity of usable galaxies from the Canada-France-Hawaii- 
Telescope (CFHT) data. We display in Figure [2]; a) the 
distribution of these source galaxies. Comparison of 
the source distribution with the exposure variation (Pig- 
ure[2](b)) hints at the correlation between the two maps. 
Nevertheless, this correlation is not to the extent that we 
worry about the impact of the number density fiuctua- 
tion on our mass reconstruction. On a ~ 1' scale (roughly 
the side of WF CCD), the peak-to- valley variation is at 
the - 7% level. 

3. WEAK-LENSING ANALYSIS 
3.1. Substructures 

We reconstruct the mass distribution using the max- 
imum entropy /maximum likelihood (MEML) algorithm 
detailed in Jee et al. (2007b), which starts with a flat 
prior and then updates the prior as the model improves. 
We performed sanity checks with three other available 
codes (Kaiser & Squires 1993; Fischer & Tyson 1997; 
Lombardi & Berlin 1999) to examine if any pronounced 
substructure is also found in the other mass reconstruc- 
tion algorithms. Figure [3] compares the results. The 
structures well inside the ooundary (i.e., the footprint 
of the WFPC2 pointings) are all visible in different al- 
gorithms whereas the structures near and outside the 
boundary depend on the reconstruction algorithms. We 
observe that the MEML reconstruction minimizes the 
boundary effects most successfully. We also claim that 
the relative amplitudes for different mass peaks are most 
reliably recovered in this reconstruction because the algo- 
rithm does not pre-smooth the galaxy ellipticity and also 
it properly distinguishes shear (7) from reduced shear 
[7/(1 — k)] through iteration (important where the pro- 
jected density is high). Therefore, hereafter our discus- 
sion on the substructure is based on the result from the 
MEML reconstruction. 

We overlay the mass contours on the pseudo-color com- 
posite of A520 in Figure |4j As in M07, we represent the 
intensity of the Chandra X-ray emission in red. The 
overall large-scale structure closely resembles the one re- 
ported in M07, although the current HST data enable 
us to resolve the substructures with much higher signifi- 
cance. We label the six mass peaks following the scheme 
of M07 and refer to them as Pl-6 hereafter. 

The most remarkable consistency between the current 
and M07 mass maps is the presence of the strong "dark 
core" (P3), which coincides with the peak of the X-ray 
emission but lacks luminous cluster galaxies. PI, P2, and 
P4 identified in M07 are in good spatial agreement with 
the equivalent peaks. 

P5 and P6 are new substructures identified in the cur- 
rent WFPC2 analysis. P5 was reported missing in M07, 
and this absence of a distinct mass peak despite the 
apparent concentration of cluster galaxies was consid- 
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Figure 1. Influence of CTI on mass reconstruction, (a) CTI pattern within WFPC2 CCDs. It shows how typical circular WFPC2 PSFs 
are distorted due to CTI. The size and orientation of the "whiskers" represent the magnitude and direction of elongation, respectively, (b) 
CTI-induced ellipticity pattern when the observational footprint is considered, (c) CTI-only mass reconstruction. The periodic variation 
due to the observational dither pattern is clear. However, even the maximum amplitude of this uncorrected systematics for the average 
source galaxy is an order of magnitude smaller than the average lensing signal. 

ered another unusual feature of M07 mass reconstruction. 
The new mass map has therefore resolved this apparent 
discrepancy. P6 coincides with some of the bright clus- 
ter galaxies, and the M07 mass map also indicated some 
overdensity in this region, albeit at a lower significance. 



3.2. Mass, Luminosity, and Significance Estimation 

To estimate the cluster mass, we combined the shape 
catalogs from the current WFPC2 and the previous 
CFHT images by using WFPC2 shapes in the overlap- 
ping region and CFHT shapes for the rest. The result- 
ing mass reconstruction is shown in Figure [5] Since the 
WFPC2 images resolves fainter galaxies at nigher red- 
shift, the amplitude of the lensing signal in the HST 
imaging is higher when no redshift correction is applied. 
Although the proximity of the cluster (z = 0.201) makes 
the amount of correction small (< 6%), we t ake this into 
account in the mass estimation (see j |3.3| for details). 
We find that the tangential shear profile at large radii 
(r > 200") is well described by a singular isothermal 
sphere (SIS) with a,, = 987 ±49 km s~^, which is consis- 
tent with the estimate of 1028±80 km s^^ from M07. We 
estimate that the aperture mass (Fahlman et al. 1994) 
within r = 710 kpc is (4.47 ± 0.48) x lO^'' Mq, again 
statistically consistent with the previous measurement 
(5.00 ±0.55) X IQi-* Mq. 

For the determination of the i3-band luminosity, we 



performed synthetic photometry using the Kinney et al. 
(1996) spectral templates, and derived a transformation 
to convert g-band (CFHT) magnitude and g~r color to 
the rest frame S-band luminosity. The total luminosity 
within the r = 710 kpc aperture is 1.7 x IO^^Lbq, which 
gives a total mass-to-light ratio 263 ± 28 Mq/Lbq- The 
updated mass, M/L, and gas fraction for the substruc- 
tures Pl-6 are listed in Table 1. 

The significance of the substructures is estimated with 
the following method. The background level cannot 
be determined reliably within the WFPC2 field (Fig- 
ure |4]) where the cluster contribution is non-negligible. 
Thus, we created a wide-field convergence map from the 
WFPC2 + CFHT data and determined the background 
level in the 13' < r < 17' annulus. The rms value of the 
convergence at the location of the peak can be estimated 
in two ways. First, we can perform a bootstrapping anal- 
ysis and compute the standard deviation with respect to 
the mean. Second, the Hessian matrix (whose elements 
are the second derivatives of the likelihood function) can 
be utilized using the Gaussian approximation of the er- 
ror distribution at the peak of the posterior distribution. 
The first method is not a viable option in the current case 
because our high-resolution MEML mass reconstruction 
requires a significant CPU time (~ 1 day with 24 cores) 
to reach the final solution. Therefore, in this paper we 
utilize the Hessian matrix to estimate the rms values of 
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source galaxy distribution exposure map 




Figure 2. Source galaxy distribution vs. exposure map. Darker shade represents higher value, (a) Source galaxy positions are smoothed 
with an FWHM=10" Gaussian kernel, (b) Shown here is the exposure map for the combined WFPC2 F814W image. The comparison of 
the two figures hints at the correlation lietween the image depth and number density of source galaxies. However, the correlation is not 
strong. On a 1' scale (the size of WF CCD), the peak-to- valley source galaxy number density variation is at the ^ 7% level with respect 
to the mean density whereas the exposure time varies significantly from mytildeA, 400 s to 13, 200 s. Therefore, the spatial variation of the 
source galaxy number density due to the exposure time variation does not cause any spurious substructures in our mass reconstruction. 




Figure 3. Mass reconstruction in A520 with different algorithms. The dashed line shows the footprint of the WFPC2 observations. The 
structures well inside the boundary are all visible in different algorithms whereas the structures near and outside the boundary depend on 
the reconstruction algorithms. 
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Figure 4. Mass reconstruction in A520. The intensity of the diffuse Chandra emission is depicted in red. The background is the pseudo- 
color composite created from the CFHT r and g passband images. The WFPC2 observation footprint is shown in orange. The white 
contours represent the convergence and the spacing is linear. The lowest contour corresponds to the ~ 2.6cr significance. The numbers 
(1-6) indicate the significant mass peaks. 




Figure 5. Large-field mass reconstruction based on the combined {HST and CFHT) catalogs. On the left-hand side, we overlay the mass 
contours on the smoothed rest-frame _B-band luminosity distribution of the cluster (linear scale). On the right-hand side, we illustrate the 
distribution of the high (red diamond, Sv^f > 1700 km s~^) and low (green diamond, Sv^f < —1500 km s~^) velocity groups. 
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Table 1 

Mass and Luminosity Properties of Substructure (r < 150 kpc) 



Substructure 


a, 5 

m s o///^ 


Aa, AS 
("/') 


Projected Mass 
(fcfoho^M©) 


Luminosity 
(fcfo'lO"iB0) 


M/L 
[HjoMq/Lbq) 


fa 


PI 


(04 54 20.76, +02 57 38.4) 


(4.4,2.7) 


2.63 + 0.48 


1.54 


171 + 31 


< 0.06 


P2 


(04 54 15.02, +02 57 09.2) 


(4.0,6.5) 


3.83 + 0.42 


3.58 


106+ 12 


< 0.08 


P3 (dark core) 


(04 54 11.07, +02 55 35.3) 


(6.7,6.5) 


4.00 + 0.38 


0.68 


588 + 56 


< 0.14 


P4 


(04 54 04.32, +02 53 51.0) 


(5.1,6.9) 


3.64 + 0.45 


2.95 


123+ 15 


< 0.08 


P5 


(04 54 16.53, +02 55 26.7) 


(6.5,6.4) 


3.03 + 0.40 


2.12 


143 + 19 


< 0.05 


P6 


(04 54 08.85, +02 53 50.2) 


(9.6,6.7) 


3.33 + 0.40 


1.23 


270 + 33 


< 0.06 



Note. — The positional uncertainty is estimated from bootstrapping. We estimate tiie aperture mass based on the method of Fahlman ct al. 
(1994). The mass uncertainties are evaluated from 1000 Monte-Carlo realizations. The gas fraction fg is derived using Caucliy-Schwartz method in 
M07 



the mass pixels. Bridle et al. (1998) demonstrated that 
the errors derived from their Hessian matrix are consis- 
tent with the values obtained from Monte Carlo simu- 
lations. We also confirmed their claim by performing 
independent simulations using a low-resolution (20 x 20) 
mass grid. We did not neglect the off-diagonal elements 
of the Hessian matrix because the mass pixels are highly 
correlated through the maximum-entropy regularization. 

There is an ambiguity in determining the absolute con- 
vergence value at the location of the mass peak because 
of the mass-sheet degeneracy. However, our wide-field 
mass reconstruction is performed in such a way that the 
convergence field in the 13' < r < 17' annulus is close 
to zero. Consequently, the convergence near the cluster 
mass peak remains virtually unchanged when we apply 
the K — >■ Ak -|- (1 — A) transformation. 

Within r < 150 kpc, the significance (defined by the 
background-subtracted convergence divided by the rms) 
of the "dark core" is ^ 12 a. This significance is slightly 
higher than the value obtained from our aperture mass 
densitometry, which gives ~ 10 cr. The difference in part 
comes from the fact that the aperture mass densitometry 
uses a limited range of tangential shears. Nevertheless, 
we find that the ratio between the aperture mass and its 
uncertainty is a conservative measure of the significance 
of the substructures (Table 1). 

Another quantity related to the significance, but a dif- 
ferent measure of the reliability of the mass peaks is their 
positional uncertainty. We performed a bootstrapping 
analysis by generating 1000 noise realizations and de- 
termine the centroid of each mass peak. Because our 
maximum-entropy method is slow and thus not practi- 
cal for this experiment, we used the Kaiser & Squires 
(1993) algorithm. The direction of the shifts during this 
bootstrapping run is nearly isotropic, and we show the Icr 
value of the distribution in Table 1. The typical la devia- 
tion is ^ 8'.' 5 corresponding to ~28 kpc at the cluster red- 
shift. Note that since the direct inversion by the Kaiser 
& Squires (1993) is somewhat noisier than the MEML 
mass map, the estimated centroid errors obtained in this 
way are likely to be overestimated (we also find that the 
significance of the dark core reduces to ~7a when we use 
these 1000 noise realizations). Nevertheless, these values 
are still small compared to the size of the mass peaks, 
which eliminates the possibility that the dark core is the 
result of any catastrophic centroid shifts of nearby peaks. 

3.3. Source Redshift Estimation 



J] 




Velocity relative to cluster mean (km/s) 

Figure 6. Galaxy redshift distribution in the A520 field. The dis- 
tribution is inconsistent vfith a single Gaussian, and the redshift 
distribution shows a possible presence of high/low velocity groups. 
However, we note that the spatial distribution of the galaxies be- 
longing to the high/low velocity groups is not concentrated near 
the dark core, but instead scattered across the cluster field (see the 
right panel of Figure p[. The shaded histogram shows the redshift 
distribution of the galaxies within the r ~ 150 kpc radius of the 
dark core. 

We select objects in the 22 < F814W^ < 27 range as 
source galaxies. For the 22 < F814W < 24 sources, 
we removed the red sequence galaxies defined by the 
CFHT data and also the spectroscopically confirmed 
cluster members. The seeing of the CFHT A520 im- 
ages is excellent, ranging from FWHM=0'.'50 to 0'.'57. 
Hence, the cross-identification between the WFPC2 and 
CHFT images can be done with high fidelity for this 
relatively bright subsample. In the magnitude range 
24 < F814:W < 27, we did not attempt to remove cluster 
contamination. At fainter magnitudes, most of the mem- 
bers will be blue, and thus the color-based cut becomes 
ineffective. Also, in this regime the blue cluster galaxy 
contamination is low because the number of background 
galaxies overwhelms that of the cluster galaxies, which 
is verified by our comparison of the number density in 
the WFPC2 field with those from the Ultra Deep Field 
(UDF; Beckwith et al. 2003) and Great Observatories 
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Origins Deep Survey (GOODS; Giavalisco et al. 2004) 
data. 

The redshift distribution of our source galaxies is esti- 
mated using the pubHcly available UDF photo-z catalog 
(Coe et al. 2006). The method is explained in detail in 
our previous publications (e.g., Jee et al. 2005a; Jee & 
Tyson 2009). The critical lensing density Scrit is pro- 
portional to (/3) — (max(0, D/s/Dg)), where D/s and Dg 
are the angular diameter distances between the lens and 
the source, and between the observer and the source, re- 
spectively. For the galaxies available within the WFPC2 
field, we obtain (/3) = 0.64. Mahdavi et al. (2007) re- 
ported that (/3) is ~ 0.60 from their analysis with the 
Ilbert et al. (2006) photometric redshift catalog. There- 
fore, when we combine the WFPC2 and CFHT shapes 
for the above mass estimation, we need to correct for this 
~6% difference in Scrit- 

4. NATURE OF THE DARK CORE 

Now that our high-resolution weak-lensing analysis 
confirms the dark core in A520, more extensive efforts 
should be made to understand the nature of this pecu- 
liar substructure. In this paper, we review and extend 
the scenarios that may lead to this significant feature. 

Compact High M/L Group. Although we cannot find 
any giant elliptical galaxies near P3, there are at least 
~ 11 spectroscopic members within a r = 150 kpc aper- 
ture. The smoothed light distribution (left panel of Fig- 
ure [5| shows a marginal indication of a faint group coin- 
cident with the X-ray peak. However, if this is a group, 
the M/L has to be extremely high, although the current 
M/L estimate (588 ± 56 Mq/Lbq) is somewhat lower 
than the previous value (721 ± 179 Mq/Lbq). One may 
suspect that the high concentration of the plasma near 
P3 somewhat contributes to this high M/L. Because 
many lines of evidence (e.g., bow-shock feature) suggest 
that the concentration of the gas originated from other 
substructures (e.g., P2 and P4), it is worth examining the 
M/L value for P3 in the absence of the gas. Subtraction 
of the X-ray gas mass from the total mass reduces the 
M/L value further to ~ 510 Mq/Lbq when the upper 
limit 0.52 x 10^^ Mq of the plasma mass in M07 is as- 
sumed. Nevertheless, this M/L value is still substantially 
higher than the mean value of rich groups. For example, 
Parker et al. (2005) quote 195 ± 29 h^MQ/LBQ for 
rich groups from their weak-lensing study of 116 groups. 
However, considering that the distribution of the M/L 
values of groups is still poorly known, we believe that it 
is premature to rule out this possibility. 

Contribution from Neighboring Substructures. M07 in- 
vestigated if a superposition of the two dark matter halos 
on P2 and P4 can lead to the detection of P3. When we 
repeat the analysis with the current updated mass of the 
two halos and the addition of P5, the resulting aperture 
mass within a r = 150 kpc would be ~ 3 x 10^^ Mq, still 
an order of magnitude smaller than what is required to 
produce the lensing signal. 

Distant Background Cluster. An extremely massive 
cluster at significantly high redshift [z > 1) can make its 
member galaxies difficult to be identified in the optical 
images, but can still signal its presence by distorting even 
higher redshift galaxies. However, the redshift depen- 
dence of the lensing signal is inconsistent with the object 
being at such high redshift. For example, when our mass 



reconstruction is performed with relatively bright source 
galaxies (thus mostly at z < 1), the mass peak is still 
visible. 

Ejection of Bright Calaxies. With A^-body simula- 
tions, Sales et al. (2007) find that a significant fraction of 
satellite galaxies are ejected during their first approach 
to the main system via three-body encounter. M07 con- 
sidered the possibility that the dark core in A520 became 
devoid of bright cluster members by the same dynami- 
cal process. This hypothesis is in part motivated by the 
presence of the rich galaxy group (P5) ^ 300 kpc east 
of the dark core. There was no significant mass associ- 
ated with P5 in M07. However, as mentioned in Sj3j our 
new analysis reveals the clear peak on top of the galaxy 
group. Moreover, the ejected component in the three- 
body encounter is always the weaker system of the pair 
infalling to the main system. Therefore, we rule out this 
possibility for the peculiar structure of A520. 

Collisional Deposition of Dark Matter. The X-ray 
analysis of the prominent bow shock feature (Marke- 
vitch et al. 2005) indicates that perhaps the subclus- 
ter (P4) was passing through the main cluster's core at 
a supersonic speed. If we attribute P3 to a deposition 
of collision-stripped dark matter, indeed A520 is a coun- 
terexample to the Bullet Cluster. Assuming that PI, P2, 
P4, and P5 equally contribute to the total mass of P3, 
M07 estimates o-^m /moM ~ 3.8 ± 1.1 cm^g"-'^, which 
exceeds the upper limit aoM /thdm < 1 cm^ g~^ de- 
rived from the Bullet Cluster (Markevitch et al. 2004), 
although we note that Williams & Saha (2011) suggest 
that a kpc-scale separation between stellar and dark mat- 
ter components in the cluster A3827 may be evidence for 
dark matter with a non-negligible self-interaction cross- 
section. Our current WFPC2-based study reduces the 
mass uncertainty of the substructures by more than a 
factor of 2 while the mass estimates remain close to the 
previous measurements. Therefore, the current improve- 
ment in precision only increases the significance of the 
above discrepancy. 

Filament along the LOS direction Numerical simula- 
tions show that galaxy clusters form where filaments in- 
tersect. Of course, the most apparent filament in A520 
might be the NE-SW structure traced by galaxies, X- 
rays, and mass. If another filament associated with A520 
is oriented along the LOS, its center should perhaps be 
near P3, because it corresponds to the approximate cen- 
ter of the large scale structure even with P3 excluded. 
As M07 noted, this filament should be sparse enough 
that no strong concentration of X-ray gas nor galax- 
ies is observable. From the velocity field analysis with 
167 spectroscopic cluster members, Girardi et al. (2008) 
report possible existence of high velocity group, whose 
rest-frame LOS velocity is ~ 2000 km s~^ relative to the 
main system, which can be viewed as indicating a LOS 
filament. 

Our analysis with 236 redshifts obtained from 
the Canadian Network for Observational Cosmology 
(CNOCl; Yee et al. 1996) redshift survey and our in- 
dependent DEIMOS redshift survey of the cluster field 
also supports the possibility of the aforementioned high 
velocity group along the LOS direction. Figure [6] shows 
the redshift distribution in the A520 field. The presence 
of the distinct bump at Vrf ^ 2000 km s~^ is clear. More- 
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over, the low end of the distribution possesses marginal 
indication of a possible presence of a low velocity group 
at Vrf ~ —1500 km s~^. 

However, the spatial distributions of both high and 
low velocity groups are not compact, but rather scat- 
tered, roughly following the large scale structure (see the 
right panel of Figure l5| , which poses a challenge to the 
scenario that these groups are responsible for the detec- 
tion of the dark core. If both groups are dynamically 
associated with the main cluster, the high and low ve- 
locity groups might be loose foreground and background 
groups, respectively, in-falling to the main system along 
the LOS direction. The velocity dispersion of the high 
velocity group is 418 ± 63 km s~^ estimated out of 34 
galaxies whose rest-frame velocities are > 1500 km s~^ 
from the center of the main system. If the group's mass 
is assumed to follow a SIS, the projected mass within an 
r = 150 kpc aperture is ~ 2 x 10^"^ Mq. Therefore, if 
we further assume that the low velocity group contains a 
similar mass, the contribution from the two LOS groups 
and the X-ray gas (< 0.5 x 10"'^^ -^0) can marginally add 
up to the required mass of the dark core. As a mat- 
ter of course, one questionable assumption is that the 
cross-sectional mass profile of the filament has a steep 
gradient, the plausibility of which requires detailed nu- 
merical studies. The galaxy distributions of the high and 
low velocity groups do not show any indication of such 
a gradient. If the hypothesized LOS filament possesses a 
smooth cylindrical mass distribution (i.e., without cen- 
tral cuspiness), it is not likely to lead to the detection 
of any isolated substructure as the dark core. Another 
immediate problem is that the filament should have in- 
creased the M/L values of other substructures as well, 
inconsistent with our measurements (Table 1). 

Aside from the above diffuse filament scenario, one can 
also imagine a thin filament coincident with P3. In this 
case we do not associate the aforementioned high veloc- 
ity group with this narrow filament. Instead, we suggest 
that the 11 cluster galaxies mentioned in the context of 
the compact high M/L group represent group-like sub- 
structure embedded in this narrow filament. This mass 
configuration can explain the lensing detection of the 
dark core, but it too poses concerns. First, while the 
additional mass in the filament would result in a high 
M/L for the substructure, an M/L of ~ 510 Mq/Lbq 
is considerably larger than the typical observed M/L of 
^ 300 Mq/Lbq for filaments (Eisenstei et al. 1997; 
Schirmer et al. 2011). Second, the hypothesis is hard to 
prove (or disprove) given the current statistics derived 
from ^11 galaxies. 

5. SUMMARY AND CONCLUSIONS 

We have presented our iJS'T/ WFPC2 lensing study 
and confirmed the reality of the peculiar mass structure 
of A520. More than a factor of three increase in the num- 
ber of usable galaxies per unit area enables us to trace 
the complex mass distribution of the cluster with higher 
resolution and with greater significance. The "dark core" 
that is coincident with the location of the X-ray luminos- 
ity peak but is largely devoid of bright cluster galaxies 
is clearly detected at the > lOcr level. With the current 
data, we reviewed several scenarios which may explain 
the presence of this unusual mass distribution. 

Dark matter self-interaction cross-section must be at 



least ~6 (7 larger than the upper limit ~lcm^g~^ deter- 
mined by the Bullet Cluster observation. Therefore, it 
is difficult to attribute the feature to dark matter self- 
interaction without falsifying the weak lensing analysis 
of the Bullet Cluster (Clowe et al. 2006), which does 
not show any significant mass clump between the two 
dominant mass peaks. Interestingly, when Bradac et al. 
(2006) combined the strong- and weak-lensing signals, 
they revealed a non-negligible mass concentration coin- 
cident with the main X-ray peak. Therefore, it might be 
worth investigating further the reality of this feature of 
the Bullet Cluster in detail. We also note that the esti- 
mate ~3.8 ± 1.1 cm^g~^ by M07 is still within the upper 
limits set by other dark matter coUisional cross-section 
studies (e.g., Natarajan et al. 2002; Hennawi & Ostriker 
2002; Bradac et al. 2008; Merten et al. 2011). 

The presence of a galaxy group with unusually high 
M/L might be a plausible solution, although the required 
M/L is considerably higher than the mean value for very 
rich groups. Since the distribution of the M/L values of 
galaxy groups is still poorly known, we cannot exclude 
this possibility yet. 

We also considered a scenario, wherein a fortuitous su- 
perposition of an LOS filament is located near the dark 
core. This hypothesis is in part supported by the spectro- 
scopic data. If the velocity dispersion of the high velocity 
group is representative of its mass and if the low veloc- 
ity group contains a similar mass, it is possible that the 
sum of the two masses can cause the lensing detection 
of the dark core. However, the spatial distributions of 
both groups are very broad and not concentrated on the 
dark core. Hence, it is difficult to rationalize that the 
projected density profile of the hypothesized diffuse LOS 
filament possesses such a cuspiness at its center. 

Alternatively, we also discussed possible existence of a 
thin filament coincident with the dark core. Although 
this might explain the lensing detection of such a pecu- 
liar substructure, the current data (only 11 spectroscopic 
redshifts) do not provide sufficient statistics to convince 
us of the presence of such a thin, but long filament near 
the dark core. More extensive spectroscopic surveys near 
the dark core area are required to test the hypothesis. 

Despite our solid confirmation on the presence of the 
dark core, we conclude that it is yet premature to single 
out the most probable cause of the dark core from the 
above scenarios. 
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